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INFLUENCE OF LEADING-EDGE SUCTION ON LIFT -DRAG BATIOS 
OF WINGS AT SUPERSONIC SPEEDS 
By Clarence B. Cohen 


SUMMARY 

A method based on linearized theory is presented for calcinat- 
ing the theoretical suction force at the subsonic leading edges of 
a f amily of wings at supersonic speeds. The method is used to 
determine the optimum sveepback angles for the tips of trapezoidal 
wings and to determine the effect of curvature of the tip contour 
on the lift-drag ratio of wing regions influenced by the tip. The 
effect of skin friction is included. The possible gain in lift- 
drag ratio from proper tip design of trapezoidal wings increases 
as the sveepback of the wing is increased. Results indicate that 
appropriately curved tip boundaries will give higher lift-drag 
ratios in the region affected by the tip than the best trapezoidal 
wing. 


INTRODUCTION 

The two-dimensional theory of thin airfoils at subsonic speeds 
gives a resultant force normal to the wing surface if only pressure 
forces on the wing surface are considered, which implies a net drag. 
For this theory, however, the velocity components became infinite 
at the leading edge. Evaluation of the effect of this singularity 
yields a suction force that exactly cancels the drag and thus sat- 
isfies the momentum requirement that the resultant force vector be 
normal to the direction of flight (reference l). 

The procedure used to derive the suction force is not 
restricted to two-dimensional wings and indicates that such a 
force may occur along any edge at which the normal component of 
the velocity becomes infinite. In the linearized .theory for thin 
wings flying at supersonic speeds, such a singularity exists along 
leading edges swept behind the Mach angle. It has therefore been 
proposed by Brown (reference 2) that a leading-edge suction force, 
which tends to counteract the wave drag, exists along subsonic 
leading edges of wings flying at supersonic speeds. The magnitude, 
of this force is calculated in reference 2 for triangular wings 


2 


MCA TN No. 1718 


swept "beyond the Mach angle. The reduction in drag predicted, "by 
the concept of leading-edge suction force has been shown "by 
Vincent! to occur experimentally for triangular wings swept beyond, 
the Mach angle. 

The perturbation-velocity components required to evaluate the 
local leading-edge suction force may be obtained by the methods of 
reference 3 for the class of wings or wing tips for which the for- 
ward Mach line from the subsonic leading edge intersects a super- 
sonic leading edge. An analysis of the suction force for such 
wings was made at the MCA Cleveland laboratory in 1948 to deter- 
mine its effect on the lift- drag ratio obtainable with various tip 
contours. 
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SIMBOLS 

The following symbols are used in this report: 
constant defining location of trailing edge of wing 
coefficient of drag due to friction 
lift coefficient 
pressure coefficient 
chord of wing 
drag force 

suction force (in flight direction) 
quantity defining cp n 

constant defining angle of straight wing edge 
lift force 

Mach number (free stream) 
distance normal to leading edge 
dynamic pressure, ^pU 2 
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▼-Tjfag) 7 

u-U2 ” u-Ug 


E distance ratio, 

S Integration area 

s distance along wing edge 

U free-stream velocity 

u, v oblique coordinates 

7 distance (in v- direction) between subsonic find supersonic 
leading edges 

x, y Cartesian coordinates 
a angle of attack, radians 

P cotangent of Mach, angle, Vm^-1 

e semiapex angle of triangular wing 

9 angle between wing edge and flight direction (positive 
counterclockwise), 

•% I 

(i. Mach angle, sin“ x ^ 

p mass density of air 

cp perturbation-velocity potential 

Subscripts: 

1 supersonic leading edge 

2 subsonic leading edge 

f friction 

n indicates direction normal to wing leading edge 
t tip point of wing (at greatest span) 

w wave 


• x, y indicates velocity components in x, y directions 
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DERIVATION OF EXPRESSION FOR LEADING-EDGE SUCTION FORCE 


In subsonic airfoil theory, the form of the equation of the 
velocity component normal to the leading edge (reference l) is 


*n 

nn»0 


G 

SE 


(i) 


where cp n = ^ is the perturbation velocity in the n direction 

(fig. 1). Thus cp n approaches infinity as n approaches zero. 

Analysis of this singularity results in a suction force per unit 
length, which can he expressed in differential form as (refer- 
ence 1) 


dF 

n 

cLs 



( 2 ) 


where s is distance along the wing edge (fig. 1). 


Equations (l) and (2) may he rewritten, as in reference 2, t d 
Include the effects of compressibility. 


Cp a 

n->0 


G 1 

s* J Zg 


(la) 


and 



= *oS 2 


(2a) 


where is the component of the Mach number nor mal to the edge. 

For subsonic two-dimensional wings, the factor G is such 
that the force defined by equation (2a) is Just sufficient to 
cancel the drag resulting from the flat-plate analysis. For 
supersonic wings, the factor G must be deter mine d from the 
appropriate value of cp^ at the subsonic leading edge. From 

reference 3, the perturbation-velocity components in the x and y 
directions near the leading edge can be written (appendix A) as 
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q> » ISt 

n->0 


m Ua 
nr->0 



> 


( 3 ) 


The coordinate system used herein is shown in figure 1. The 
relations "between the oblique and Cartesian coordinates are given 
by 


u = H (x-Py) 


x = | (v-rn) 


v = ~ (x+py) 

y - 5 (t “ u) 


(4) 


Equations (3) may be combined to obtain the normal velocity 
component near the edge 


» Hoc, 
n->0 


where 





E = 


r-VjCug) y 
U-Ug U-Ug 


(5> 


The geometrical relation between E and n as n-»0 is 
shown in figure 2, where 

n = PQ sin (0-Hi) = — sin (0+p.) 


( 6 ) 


■ following expressions are derived from the relations 
between the oblique and Cartesian coordinates: 
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With the aid of equations (6) and (7), equations (la) and (5) 
may he solved for G 2 . 



This value of G substituted into equation (2a) yields the 
normal force per unit length along the wing edge. 



O) 


With the use of the relation 


dF ds 
ds dv 
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there results 



The leading-edge suction force can now be expressed in Inte- 
gral form as 



The quantity (-5)4* in equation (ll) is given in 

figure 3 as a function of 6 for three values of p. With the 
aid of curves such as these and. a drawing of the wing in question, 
equation (ll) may be evaluated by a numerical summation. 

The force defined by equation (ll) is obtained when the radius 
of curvature of the leading edge is vanishingly small, and when the 
velocity components normal to the edge become infinite at the edge. 
In practice, however, the Infinite velocity components predicted by 
linearized theory correspond to a high negative pressure coefficient 
near the edge. By rounding the subsonic leading edge to increase 
the projected area on which the lower pressure can act, it is 
possible to reduce the total drag of the wing. As the edge radius 
is increased, an optimum curvature that corresponds to a maximiza- 
tion of the product of area times suction force per unit area, 
should be obtained. It has been shown experimentally that round- 
ing the leading edge of triangular wings swept beyond the Mach 
angle increases the lift-drag ratio, but that continued rounding 
beyond a certain value produces no further gains. 


APPLICATION OF LEADING-EDGE SUCTION-FORCE EQUATION 


For a given main surface (fig. 1), the optimum wing is one 
whose tip has the greatest lift-drag ratio L/D. If the tip has 

1 

Xa Xi CD 

no subsonic edge, this ratio is simply - 


D By. + Dp 


1 + 


Dp 


where 
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is the -wave drag the skin-friction drag. If the tip 

Tift« a subsonic leading edge its lift-drag ratio becomes 


L L 

D D + D ~ - F 
w x 



+ % _ F_ 
•^w 


( 12 ) 


where F is the suction force at the subsonic leadi n g edge. The 


tip with the minimum value of 


(r-4) 


therefore has the highest 


lift-drag ratio. If frictionless flow is assumed, the param- 
eter F/H w should be maximized. 

As an example of the procedure for obtaining wings of optimum 


plan boundaries, the quantity 



w ill be derived for a 


wing with a straight-line supersonic leading edge v^= k-jU. The 
quantity V in equation (ll) then becomes v - k-jUj (fig. 4) and 
the suction force is 



where v t is the value of v at idle extreme tip and u = Ug(v) 
is the equation of the tip edge. 

The wave drag for a wing tip of this type is (appendix B) 




\ 

< 

i 


and the friction drag is 


% = \ X (tip area) = 


nn_ .«rr2 

— V,X^'’ 


r* 


§-< T - k i u 2 )_ 


dv 


(15) 


mV%A AAMA^nWkJl V 4«> - 


u<%uu u xxi ww o^ucioxLwx lor uuo LTA xiing eage aexermiiisa Tine onora or me x n-ng 
(fig. 4). Then combining equations (IS) to (15) 



If the wing-tip equation is kept in the general form ug(v), the conditions for min- 
imizing equation (16) can only be indicated. The procedure is to set the differential of 
equation (16) equal to zero; the solution to the resulting integral equation^ however; 
not yet been determined. 


If the wing is assumed to have a straight tip edge given by u^, = kgV, then the plan 

form 1b a trapezoid and the integrals in equation (16) can be evaluated (appendix C) to - 
yield 
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where 


\ “ (3% - atjkg - 1) 

- ® < k.. <0 

0 < Itg < 1 


(17) 


and k^ and kg determine the sweep angles of the leading and 
trailing edges 0^ and 6 g, respectively. The angles corre- 
sponding to several values of the constants are: 


k l 

0 1 

k 2 

e 2 

0 

-v- 

0 


-1 

«/2 

1 

0 

-09 

Ji 




The equations relating the angles and the constants are 



CfcL - 1) 1 
»(*! + 1) 


0 


2 



(1 - kg) 

P(1 + kg) 


> 


, (18) 


As an indication of the variation of leading-edge suction 
force for this type of wing, the effect of 0g upon f/D^., as 

calculated from equations (17) and (18) with C D f = 0, is shown 
in figure 5 for 0-^ = 65° (k^ = -2.75) at M = fj2. An optimum 

value of 0g occurs at about 20°. 

In order to obtain an expression for optimum values of 0g 
(or kg) as a function of k-j_, equation (17) must be differ- 
entiated and set equal to zero. When differentiated with respect 
to kg, the equality yields 
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2kg 2 k x + kg (3-7k 1 ) + 3k 1 - 1 + = o (19) 

or ^ 

Differentiation with, respect to k^ yields no true rnnytimim for the 
type of wing treated. The variation of the optimum value of kg 
with ki, as expressed by equation (19), is shown in figure 6 for 

PC D f 

several values of — . 

or 

/ 

A wing with 0^ = 65° and. three alternate tip curves A, B, 

and C passing through points 0 and t is shown in figure 7. Tip 
curve B is swept so that 0 g = 20° (the optimum for this trape- 
zoidal wing from fig. 5). Curves A and C are arbitrary curves, 
symmetrical about B, drawn through 0 and t. 

The following values of F/D^ were obtained by numerical 

integration of equations (13) and (14) for these tips; the result- 
ing effect on lift-drag parameter aL/D for frictionless flow is 
also shown: 


Wing 


oL/D 

(frictionless 

flow) 

A 

0.39 

1.64 

B 

.31 

1.45 

C 

.20 

1.25 


These results indicate that the lift-drag ratio may be Increased by 
the use of appropriately curved wing boundaries at supersonic 
speeds. 

As an indication of the magnitudes of F/Dy theoretically 
obtainable with other type plan forms, the variation of f/d w 

with apex angle for a triangular wing is shown in figure 8. tmh 
relation is based on the equations in reference 2. 

The variation of maximum lift- drag ratio with kj_ is shown 

in figure 9 for trapezoidal tips. These lift-drag ratios were * 
calculated from the optimum values of kg shown in figure 6. 
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The lift-drag ratio for the main surface of a trapezoidal wing is 
also shown for comparison. These curves 'show that the maximum 
lift-drag ratio of sweptback trapezoidal wings is obtained, not 
when the wing is cut off at the Inner Mach line (main- surface 
curves), but rather when the optimum tip angle, based on considera- 
tion of leading-edge suction, is chosen. (This result held for all 


values of the friction parameter 


PCp.f 


shown. ) For sweptf orward 


a. 


trapezoidal wings, however, figure 9 indicates that the optimum 
lift-drag ratio L/D is obtained when the tip is cut back at the 
inner Mach angle. 


The lift-drag ratio is shown in figure 10 as a function of 
angle of attack for the wing tips previously mentioned, for the 
main surface, and for a tip cut off in the flow direction (tip D). 
The calculations were made for M = */2 and = 0.005. This 

figure indicates that appropriately curved plan forms give more 
favorable lift-drag ratios than the optimum trapezoidal wing. For 
wings of low aspect ratio, the effect of the tip design becomes 
relatively more important; a complete study must therefore include 
the effect of aspect ratio. From the standpoint of high lift-drag 
ratio, a tip such an, D (fig. 10) is poor for any aspect ratio 
inasmuch as a higher lift-drag rat ip can be obtained merely by 
cutting off the wing at the Mach line. 


CONCLUSIONS 

The following conclusions result from a method for evaluation 
of the theoretical leading-edge suction force on a type of wing 
for which the forward Mach lines from the subsonic leading edge 
intersect a supersonic leading edge: 

1. The optimum lift-drag ratio of sweptback trapezoidal wings 
is obtained when the optimum tip angle, based on consideration of 
leading-edge suction, is chosen. For sweptf orward trapezoidal 
wings, however, the optimum lift-drag ratio is obtained when the 
tip is cut back at the Mach angle. The possible gain in lift-drag 
ratio from proper tip design of trapezoidal wings increases as the 
sweepback of the wing is increased. 

2. Higher lift-drag ratios may be obtained with wing tips 
whose edges are appropriately curved than may be obtained with 
straight-line wing tips. 


Lewis Slight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
CleveH and,, Ohio, July 15, 1948. 
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APFEHDIX A 


PEBTlIRBATION-mDCITY COMPGHEHTS 
NEAP SUBSONIC LEADING EDGE 


Sep 


The x- component of the perturbation velocity cp = where 

x ox 

<P is the perturbation-velocity potential, is given in reference 3 


as 


Sep _ Ua 
Sx 2J3jr 



dv'-du' 


/y/Xu-u’ ) (v-v‘ ) 



(-a) 


(>•» 


dv' 


a/(u-u’ ) (v-v* ) 


(Al) 


vhere the primes denote the variables of integration. 



The integrals are evaluated along the supersonic leading edge 
from A to B and along the Mach line from B to Q. When point P 
approaches the subsonic leading edge (5^ (n-*0), points A and B 

' converge and equation (Al) reduces to 
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(A2) 


which, upon integration becomes the first of equations (3) . The 
expression for cp is obtained similarly. 
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APPENDIX B 
DRAG FORCE 

In accordance with 'linearized theory, the Wave drag of a thin 
flat plate is 


Dy. = aL = aC^qS 


(Bl) 


where the lift coefficient can he expressed as 


h .iIL 


C_dS 


II 


(B2) 


as 


For a wing of the type shown in figure 4, reference 3 shows 
that the pressure coefficient can he expressed as 


2a 

itp 


(>• 


anjA Iv-kjyg l-k x .-L Mug-u) 


dvy /y u-ug ^ 


tan 




(B3) 


A combination of equations (Bl) to (B3) yields 


E 



dv 


Ja gO 


(v) 


LS) 

Hws 

1-ki 

fe(u g-u) 

[ 

1 u_u 2 

T OqJl r A 

V-kiUg 


du 


(B4) 


which upon integration with respect to u becomes equation (14). 
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APPENDIX C 


SUCTIOH FORCE AND DRAG FOR A TRAPEZOIDAL WHIG 
The equation of the tip of a trapezoidal wing is 

Ugfr) = kgV 


(Cl) 


By the use of this equation, if the value of is known, the 

wave drag, friction drag, and suction force can he obtained, from 
equations (13) to (15). For the trapezoidal wing 


T * = (i-vJ 


M 


^(k-L + l) 2 - 4k x 


23 


(C2) 


The use of equations (Cl) and. (C2) in equations (13) to (15) yields 


F - ‘§r i ?' k i * l > 2 - tt i] (ra) 


_ ouVc 2 JR(% + l) 2 - ^ll (3*1 " 2^ - l)( 

8 P 2 


(l-kjkg) k lN /=^ 


(C4) 


pU^C^j. [m 2 ^ + l) 2 - 4ki 
= 83 ^(kikg-i) 


(C5) 


By subtracting equation (C3) from (C5 ) and by dividing the result 
by ‘equation (C4), equation (17) is obtained. 
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Figure 1. - Coordinate system used in analysis of wing tip. 
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Figure 5. - Variation of F/l^ for tip of trapezoidal wing with 

0 1 - 65° at Ms ^ 2 . 
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Figure 8. - Variation of F/D w with .apex angle for triangular 
wings. (Based on equations from reference 2.) 



Maximum lift-drag parameter* alyto 



Figure 9* - Maximum lift-drag ratio for regions of trapezoidal wings When effeots 
of leading-edge suction and friction drag are considered* 
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Lift -drag ratio 
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